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ABSTRACT 

We present near-infrared emission line images of the circumnuclear ring in NGC 1068. 
We have measured the Br7 fluxes in a number of star forming complexes and derived 
extinctions for each of these by comparison with Ha. We investigate the star forming 
histories of these regions and find that a short burst of star formation occured co-evally 
throughout the ring within the last 30-40 Myr, and perhaps as recently as 4-7 Myr 
ago. The 1-0 S(l) flux and S(l)/Br7 ratios indicate that as well as fluorescence, shock 
excited H2 emission contributes to the total flux. There is excess H2 flux to the North- 
West where the ionisation cone crosses the ring, and we have shown it is possible that 
the non-stellar continuum from the Seyfert nucleus which produces the high excitation 
lines could also be causing fluorescence at the edges of molecular clouds in the ring. 
The nuclear 1-0 S(l) is more extended than previously realised but only along the bar's 
major axis, and we consider mechanisms for its excitation. 

Key words: galaxies: individual: NGC 1068 - galaxies: Seyfert - galaxies: starburst 
- infrared: galaxies. 



1 INTRODUCTION 

NGC 1068 is considered an archetypal Seyfert 2 galaxy, and 
being nearby and bright (14.4 Mpc and Lir = 2 x 10 11 L Q , 
Ringberg Standard) has been observed extensively. In rc- 
' cent years much effort has been devoted to studying the 
nucleus at ever higher spatial resolutions, exemplified by a 
direct 8.4 GHz image of the inner torus at a resolution of 

2 milli-arcsec (Gallimore, Baum & O'Dea, 1997), and adap- 
tive optics imaging (Thatte et al. 1997) which showed that 
94 per cent of the K-band light in the central 1 arcsec orig- 
inates from a ^ 30 milli-arcsec source, the remainder being 
due to a nuclear star cluster. 

But the circumnuclear ring at a radius of 15-16 arcsec 
should not be overlooked as here the Hn regions, equivalent 
to a string of M 82-type galaxies, rival the nucleus in terms 
of bolometric luminosity. Telesco & Decher (1988) explained 
the ring in terms of inner Lindblad resonances associated 
with a barred potential, and Scoville et al. (1988) showed 
that there was also a bar extending from this ring closer in 
towards the nucleus. Such a hierarchical structure is one of 
the most efficient ways in which to transport gas over the 3- 
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4 orders of magnitude in scale from the disk, not only into 
the circumnuclear regions, but into the nucleus itself. Un- 
derstanding what processes are operating in the ring, and 
how and why they do so is a necessary step towards a global 
picture of NGC 1068. Atherton, Reay & Taylor (1985) pre- 
sented a velocity resolved Hcv image of the ring and circum- 
nuclear environment, finding that there were a large number 
of clouds with highly turbulent motions, and possibly a large 
scale expansion of the disc. This expansion was also seen by 
Planesas, Scoville & Myers (1991) who made high resolu- 
tion CO observations. The CO observations revealed a large 
number of cloud complexes in the ring and nucleus, which 
appear to be much warmer (~50K) and more massive (up 
to several xlO 8 Mq) than molecular clouds in the Galaxy. 

The inner disk has one of the highest blue surface 
brightnesses known (Keel & Weedman, 1978) and is the ori- 
gin for about 1/3 as much Ha again as the ring. Bland- 
Hawthorn, Sokolowski & Cecil (1991) argue that this diffuse 
ionised medium, characterised by a high [Nii]/Hq? ratio, is 
excited by scattered nuclear radiation, whereas the high ex- 
citation lines (Evans & Dopita 1986; Bergeron 1989) arise 
in gas which has a direct line of sight to the Seyfert nu- 
cleus. The morphology of the high excitation gas has been 
studied in detail by Pogge (1988) and Unger et al. (1992), 
who traced it out beyond the ring at PA 11-51°. An ioni- 
sation cone such as this would naturally occur as a direct 
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result of collimation by a molecular torus, but simple inter- 
pretation of the observations in terms of this model is not 
possible. Unger et al. demonstrate that a proper analysis re- 
quires careful study of the complex geometries of the radio 
jet, [O ill] cone, UV field, etc. 

Because of its brightness, proximity and favourable ori- 
entation, NGC 1068 provides an excellent laboratory for 
studying the Seyfert nucleus, inner disk and circumnuclear 
ring. Such work in many different wavebands (see the col- 
lage in Bland-Hawthorn et al. 1997) is already providing 
insights into how these phenomena are inter- linked, a topic 
which is receiving considerable attention at present. One as- 
pect which has been previously omitted is the circumnuclear 
near-infrared line emission and in this paper we present the 
first such 1-0 S(l) and Br7 images, which have been made 
possible because of the wide field available with the Fabry- 
Perot and IRCAM3 on UKIRT. We study individual star 
forming complexes, and by comparison of the recombination 
line strengths determine the extinctions to each region. We 
look at the timescales for the star formation, determining 
whether it is co-eval throughout the ring and if there have 
been previous episodes of activity in Section |E| In Section ^ 
we consider triggering mechanisms with reference the effect 
of the ionisation cone and radio jet, and in Section fj] we 
discuss the origin of the 1-0 S(l) emission in the ring. Our 
signal-to-noise is sufficient to have revealed H2 emission near 
the nucleus which is more extended than was previously re- 
alised, and in Section ^| we consider what may be the cause 
of this. We summarise our conclusions in Section 9l 



2 OBSERVATIONS AND DATA REDUCTION 

Images of the 2 pirn continuum, and H2l-0S(1) and Br7 
emission lines were obtained with the UKIRT 3.8-m tele- 
scope on the nights of 19 and 20 December 1996, using the 
tip-tilt facility. A Fabry-Perot etalon, with spectral resolu- 
tion of 325 kms -1 , was used in conjunction with IRCAM3 
at the Cassegrain focus. IRCAM3 is a near-infrared camera 
housing a 256 x 256 InSb array with a scale of 0.3 arcsec 
per pixel, giving an unvignetted field of view of more than 
60 arcsec. Cooled narrow band filters (FWHM 2.4 per cent) 
were used for order sorting, and wavelength calibration was 
achieved by peaking up on the 2.1171 /im line from a Kryp- 
ton lamp. 

The observations were carried out as follows. For each 
line, 2 or 3 minute observations of the object were made 
at the on-line wavelength and also of the continuum at 
slightly shorter and longer wavelengths (shifted with re- 
spect to the line centre by ±1200 kms -1 for 1-0 S(l), and by 
-1200, +1080 kms" 1 for Br7). Sky frames were also made 
at each wavelength. This was repeated until sufficient on-line 
integration had been achieved. 

IRAF was used for all the image processing and anal- 
ysis as described below. A flatficld was made from the 
bias-subtracted 3-minute sky-frames. This was then divided 
into the object frames after each corresponding sky frame 
had been subtracted. Frames were registered on the bright 
peak of emission from the Seyfert nucleus. Each frame was 
then flux calibrated. Non-photometric conditions on the first 
night and saturation of the star used for Br7, meant that 
only the 1-0 S(l) line on the second night could be calibrated 



directly from a standard star (HR1091, an Al star with no 
spectral features at the wavelengths used). All the other 
frames were self-calibrated; that is, the continuum (sampled 
only in regions with no line emission) was scaled to match 
the calibrated frames. This was feasible because the con- 
tinuum emission in NGC 1068 is bright over a sufficiently 
extended area. 

Multiple ghost images were a severe problem with this 
object. However, as the location of the ghosts was different 
for the two nights on which the observations were made, 
these could be masked out. This was carried out while co- 
adding the individual calibrated frames. Finally, the con- 
tinuum frame was subtracted from the on-line frames to 
leave ghost-free line-only images. The total on-line integra- 
tion times were 42 minutes for Br7 and 30 minutes for 1- 
0S(1). 

The resolution (including the effects of seeing and reg- 
istering, as well as the enhancement afforded by the tip- 
tilt) was measured from the standard star profiles to be 
1.0 arcsec. Unfortunately, as the emission is extremely faint, 
the images in Figure hi have had to be smoothed with a 
median filter to enhance the signal-to-noise, reducing the 
resolution to either 1.4 or 1.7 arcsec as shown. 



2.1 Calibration & the Velocity Field 

One important point to consider is whether the calibration 
is significantly affected by the velocity distribution of the 
line emitting gas itself, since Atherton et al. (1985) showed 
that the line-of-sight velocities of Ha in the circumnuclear 
ring varied over 1020-1260 kms' 1 . For our observations we 
centred the FP at 1136 kms" 1 , corresponding to a helio- 
centric velocity of 1114 kms" 1 , so that the 325 kms -1 res- 
olution of the FP covered the range 950-1280 kms" 1 . How- 
ever, there is a shift in transmitted wavelength (towards 
shorter wavelengths) which increases to larger radius, reach- 
ing SX = 160 kms" 1 at r = 30 arcsec. Since SX oc r 2 , we es- 
timate that in the ring SX = 40 kms -1 . We have taken this 
into account when correcting fluxes for the variations in ve- 
locity. The correction required has been calculated as 1/A(v) 
where A(v) is the value of the Airy function at a velocity 
v offset from peak transmission. The velocity dispersion of 
the lines has not been included since it would require de- 
tailed knowledge of the line profiles across the entire field, 
and would only entail an additional correction varying from 
almost zero to a maximum of 10 per cent in the South- West 
where a > 135 kms -1 . 

The approximate velocities of the Br7 components 
marked in Fig. [j] are given in Table |l| together with the 
correction factor. Throughout the rest of the paper, the ob- 
served fluxes given have been velocity corrected. The whole 
image has not been corrected as, instead of simplifying its 
interpretation, it would serve only to introduce a variable 
detection limit across the field. 



3 NEAR-INFRARED LINE MORPHOLOGY 

Images of the Br7 and 1-0 S(l) line emission, and 2 /im con- 
tinuum are presented in Figure |l|. The central regions have 
not been contoured; instead the location of the nucleus, as 
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Figure 1. Contour plots of a 72 X 72 arcsec 2 field at the centre of NGC 1068. The calibration shown in the line images has not been 
corrected for the velocity field. Upper Left — Br7 line emission, with contours at 2.9, 5.7 & 8.6 xlO -19 W m 2 arcsec 2 ; the lowest 
contour is 2.5<r. Upper right - 1-0 S(l) line emission, with contours at the same intervals; the lowest contour is 3<r due to the greater 
smoothing used. Lower — 2 /im continuum after subtraction of best-fit ellipses, with contours from 0.1-2.3 mjy arcsec -2 at intervals of 
0.2 mjy arcsec -2 ; the lowest contour is 2.5(X. The central (peak) region has not been contoured and instead the location of the nucleus 
is marked by a large dot. 



defined by the centroid of the continuum emission, has been 
marked with a large dot. 

3.1 2 /j,m continuum 

Calibration of the continuum in several annuli compares 
favourably to that in Scoville et al. (1988). The emission is 
dominated by elliptical isophotes, and in order to better re- 
cover the fainter underlying structure, these have been mod- 



elled and subtracted using tasks in the ISOPHOTE package 
within IRAF. In order to preserve this structure, its inclu- 
sion in the fitting process was avoided by using a clipping 
algorithm which removed the 40 per cent of pixels with the 
highest counts before fitting each isophote. The bar is very 
clear lying along position angle 46°, as has been measured 
by other authors (eg Scoville et al. 1988). At the regions of 
peak intensity in the residual image, the bar has a surface 
brightness approximately 75 per cent of the underlying ellip- 
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Table 1. Flux Correction Factors 



detected at a much higher signal-to-noise, is consistent with 
our measurement. 



Region 


Position ( 


arcsec 


Velocity 


Correction 




from nucleus) 


kms 


factor 


A 


11 4 E 


10 2 N 


1050 


1.01 


B 


9 1 E 


12 9 N 


1090 


1.02 


C 


ft n T? 
D.U Pj 


iq ft TS.T 


1120 


1.11 


D 


3.7E 


15. 3N 


1130 


1.16 


E 


5.2W 


15. 3N 


1220 


1.90 


F 


13.2W 


7.2S 


1270 


2.57 


G 


10.6 W 


10.7S 


1260 


2.41 


H 


7.4 W 


13. IS 


1200 


1.68 


I 


4.5 W 


14.3 S 


1180 


1.49 


J 


1.7W 


13.3S 


1160 


1.33 


K 


1.7E 


13. OS 


1140 


1.21 


L 


6.8 E 


11.6S 


1060 


1.00 


M 


19.2 E 


3.4 N 


1010 


1.12 


N 


20.5 E 


9.0 N 


1030 


1.05 





20.8 E 


12. IN 


1030 


1.05 



tical isophotes. Closer to the nucleus, the clipping algorithm 
may not have been efficient in excluding pixels in the bar 
due to the small size of the ellipses and width of the bar. 
We therefore do not consider the apparent lack of structure 
close to the nucleus to be significant; indeed the H2 emis- 
sion suggests that the bar does continue much further in (see 
also Section |^) . All the remaining structure in the residual 
image, particularly at the southern end of the bar and the 
spiral arm west of the northern end of the bar, also appears 
in the original continuum image. 

It is interesting that the strongest continuum in the ring 
follows the 10 fim morphology rather than that indicated by 
Ha or Br7 knots, suggesting that the emission mechanisms 
may be associated. 



3.2 Molecular Hydrogen 1-0 S(l) 

The 1-0 S(l) emission is relatively strong at the end of the 
bar to the North-East, as well as projecting from the nu- 
cleus along the first few arcsecs of the bar. The morphology 
in the central few arcsecs is consistent with that observed by 
Rotaciuc et al. (1991) and Blietz et al. (1994). The nuclear 
region shows a strong peak almost due East of the nucleus 
and a weaker one to the South- West (Fig. Q), but these au- 
thors did not observe the more extended nuclear emission 
as their data went less deep. As a rough consistency check 
of the total flux in a 6 x 6 arcsec 2 aperture, we first interpo- 
lated across the pixels badly affected by the saturation, and 
then applied a 20 per cent correction for the 200 km s~ line 
width. The resulting flux of 13.8 ±0.1 x 10" 17 Wm" 2 agrees 
well with that measured by Oliva & Moorwood (1990) in a 
similar aperture, and also Blitz et al. (1994). 

Reay, Atherton, and Walton (1986) previously detected 
1-0 S(l) at two points in the ring using a spectrometer with 
a 5.4 arcsec aperture. At the North-East end of the bar they 
measured an intensity almost an order of magnitude greater 
than ourselves. We believe that large error bars (points on 
the spectrum are 2-3<r above their continuum level) com- 
bined with an underestimate of the continuum level, have 
led to an overestimate of the line flux. Their nuclear flux, 



3.3 Hydrogen Recombination Br7 

The Br7 clearly delineates sections of the circurnnuclear 
ring, showing distinct similarities to that seen in the Ha line 
map (Bland-Hawthorn, Sokolowski & Cecil, 1991). Although 
the strongest emission appears to be to the North-East, this 
is partly due to the effects of the velocity field as described 
in Section 2.1 . Table ^ shows that the emission to the South- 
West is as strong as that in the North-East section of the 
ring. The relative apparent weakness of the emission to the 
West is also due to velocity effects, with the result that any 
Br7 emission is below our detection threshold. In the nu- 
cleus there is evidence for a single strong peak, as observed 
by Rotaciuc et al. (1991). 



4 EXTINCTION 

The extinction towards individual Hn regions may not be 
uniform. Young, Kleinmann & Allen (1988) mapped the ex- 
tinction using the ratio of [Sm]A9532A to (Ha + [Nil]). 
They found that it peaked at A\r = 1.5 mags to the North- 
East, and at Ay = 2.5 mags to the South- West. However, the 
assumption made in order to produce these maps has been 
questioned by Kennicutt & Pogge (1990) who found an order 
of magnitude scatter in the intrinsic [S ill] /Ha ratio in Hn 
regions. They suggested that variations in the observed ratio 
might reflect variations in the level of gas excitation rather 
than dust extinction, and in the circurnnuclear regions of 
NGC 1068 several different line excitation processes may be 
operating, including photoionisation from the nucleus and 
possibly relativistic electrons associated with the nuclear jet 
(Cecil, Bland & Tully 1990). 

We have therefore measured the extinction in differ- 
ent regions by comparing our Br7 image to an Ha image 
kindly provided by J. Bland-Hawthorn (Bland-Hawthorn et 
al. 1991). This image shows pure Ha emission, as the con- 
tribution from [Nil] has been mapped and subtracted by 
fitting spectra at each pixel. The image was calibrated by 
fixing the total line flux to be that determined previously by 
these authors, 4.6 x 10 34 W m" 2 . This was further verified by 
determing the difference in (Ha+[N 11]) flux (taking the aver- 
age ratio Ha/(Ha+[N 11]) = 0.5) between a 13.5 arcsec aper- 
ture (McQuade, Calzetti & Kinney, 1995) and a 64 arcsec 
aperture (Young et al. 1988). These agreed to within 9 per 
cent. The Br7 image was smoothed to a similar resolution 
(approximately 3 arcsec FWHM) and re- binned to a scale of 
0.5 arcsec pixel -1 to match the Ha image. The images were 
then registered with reference to several prominent Hn re- 
gions in the ring. Br7 and Ha fluxes were measured in the 
regions A-O using 3 arcsec apertures. The Br7 fluxes were 
then velocity corrected as described in Section 2.1 before 



deriving the extinction. This was done assuming an intrin- 
sic ratio Ha/Br7 = 103.6 (Osterbrock 1989) and using the 
Howarth (1983) galactic extinction curve. The results are 
given in Table |2[ 

Typical values are Ay = 1-3 mag similar to those found 
by Young et al. (1988), suggesting that reddening correc- 
tions to infrared line fluxes are relatively small. This result 
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is important because by using a diagnostic in the K-band we 
should be able to detect H II regions to optical depths much 
greater than an Av of 3 mag since T2.iy,m = 1 at Av— 9 mag 
The issue cannot be avoided by considering, instead of the 
usual 'screen' model, a 'mixed' model which allows one to 
infer significantly higher Av but for which we still derive 
only Ay = 3-17 mag. It implies that in regions where we see 
star formation, either it is occuring only near the surface, or 
the molecular clouds are not highly optically thick. The lat- 
ter option seems unlikely from CO observations which imply 
considerably higher optical thicknesses for the clouds, and 
the former appears to require some external influence on the 
star formation. 



5 TIMESCALES FOR STAR FORMATION 

The spatial resolution in Fig. [j] has allowed us to measure 
the FWHM of some Br7 knots as 1.3arcsec (90 pc), and 
in the unsmoothed image we have measured the brightest 
knot to have a FWHM of 0.9arcsec (60 pc). These upper 
limits to the cluster sizes, and the masses we derive in this 
section, suggest that each region A-O, if not an individual 
OB association, is probably a relatively small number of such 
associations. 

We can determine the age of the star formation by mea- 
suring the equivalent width (EW) of the Br7 line in these 
clusters. In order to estimate the continuum flux density as- 
sociated only with young stars we have subtracted the un- 
derlying isophotes (Fig. nl). Careful comparison of the images 
suggests the labelled regions even close to the bar will suffer 
little or no contamination from the continuum in the bar. 
The Br7 and continuum fluxes in 3 arcsec apertures at these 
sites are given in Table |^, and show that EW(Br7)=20- 
200A. We use the evolutionary synthesis models of Leitherer 

6 Heckman (1995) to determine the starburst history, mak- 
ing the assumptions of solar metallicity and a standard IMF 
(Salpeter in the range 1-100 Mq). 

If the star formation is instantaneous (i.e., the duration 
is small compared to the time elapsed since it occured) then 
around the whole ring, it occured from 4-7 Myr ago. The 
similarity in the ages even with an order of magnitude dif- 
ference in EW is due to the very rapid decrease in EW with 
time, reflecting the speed of evolution of the most massive 
stars away from the main sequence. The ionising flux JVtjv 
is derived from the dereddened Br7 flux making the usual 
assumption that 70 Lyman continuum photons are required 
for each Br7 photon produced. This and the age can then be 
used to determine the mass of stars formed in each region, 
which is 10 5 -10 6 M Q . 

Alternatively, if the star formation is still progressing, 
then the ages range from 6 Myr to more than 300 Myr (un- 
less the IMF is truncated). Although the total mass of stars 
formed would be of the order of 10 6 -10 7 M© in each complex, 
this is still several times smaller than the available gas mass 
(estimated from the CO luminosity, Planesas et al. 1991) and 
so does not rule out the scenario. We believe that continuous 
star formation is unlikely as it would need to be sustained 
over timescales of more than 100 Myr in highly specific re- 
gions. In a single star cluster, once massive young stars have 
formed, the turbulence and heating from these would be- 
gin to disrupt the remainder of the giant molecular cloud 



(CMC) and inhibit further star formation. After a few mil- 
lion years these stars would evolve to produce supernovae, 
injecting more kinetic energy into the local ISM. Further 
star formation would be halted after only a few generations 
of OB stars. Blitz & Shu (1980) and Blitz (1991) used obser- 
vational and theoretical evidence to argue that GMCs which 
generate OB associations cannot have lifetimes much greater 
than 20-40 Myr. A simple persuasive line of reasoning con- 
cerns the gravitational binding energy of a CMC, typically 

10 50 erg for a 2 x 10 5 M Q cloud of size 40 pc. By estimat- 
ing the luminosity of OB stars formed over 10 7 yrs and the 
efficiency with which this could be converted into kinetic en- 
ergy Blitz & Shu found that 10 51 erg were available over this 
time, which could produce expansion speeds of ~ lOkms -1 
capable of dispersing the cloud. Thus even without consid- 
ering the kinetic energy released by supernovae (typically 

10 51 erg for a type II supernova) , the active star formation 
episode in GMCs is necessarily short lived. More detailed 
calculations of the mechanical energy deposited in the ISM 
are included in Leitherer & Heckman's models and we apply 
these to NGC 1068. They show that for a cluster undergo- 
ing constant star formation with the typical Br7 luminosity 
observed, it exceeds the binding energy of ~ 10 55 erg for 
the larger molecular clouds (Planesas et al. 1991) after 30- 
40 Myr. This therefore provides an upper limit to the du- 
ration of the star forming episodes here. The masses of the 
clouds may actually be much less as there is some evidence 
that the CO-to-Eh conversion factor has been overestimated, 
in which case the the binding energy would be significantly 
smaller and the limit on the duration of star formation could 
be less than 10-20 Myr. 

We conclude that the active phase in the ring is short- 
lived, these are young star clusters, and that the star forma- 
tion is co-eval on a timescale of 10 Myr. The 2 /im continuum 
in these regions can be accounted for entirely by this star 
formation, so that we find no evidence for previous episodes 
of star formation. Since strong Ha is observed around the 
whole ring, the arguments above can be extended to the re- 
gions where we could not detect Br7. Even to the West 
where we observe continuum emission (0.1 mjy arcsec -2 ) 
but no Br7 (< 0.22 x 10 -18 Wm -2 arcsec -2 averaged over a 
3 arcsec aperture and velocity corrected), the upper limit for 
EW(Br7) is 35A consistent with an age of 6 Myr. We there- 
fore also conclude that this is probably the first episode of 
star formation in the ring. 



6 TRIGGERING STAR FORMATION 

The existence of the circumnuclear ring in NGC 1068 has 
been attributed to gas settling between the inner Lindblad 
resonances (ILRs) as a result of the action of a barred grav- 
itational potential (Telesco & Decher 1988). In such situa- 
tions we expect vigorous star formation to occur as a direct 
result of the increased cloud density at the ILR. The Schmidt 
Law describes the star formation rate as dependent on the 
local gas surface density (SFR cx Pg as with 1 < n < 2, 
eg Kennicutt 1989) because the higher cloud density leads 
to an increased cloud collision rate so that GMCs are more 
likely to form and compression induced cloud collapse to 
occur, both of which will enhance the star formation rate. 
This appears to be substantiated, at least qualitatively, by 
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Table 2. Star Formation Properties 



Region Ay Observed flux" Corrected flux'' 





mag 


io- 18 


Wm- 2 


io- 18 


Win" 2 






Br7 


1-0 S(l) 


Br7 


1-0S(1) 


A 


1.2 


3.3 ±0.2 


3.5 ± 0.3 


3.7 ±0.2 


3.9 ±0.3 


B 


2.0 


4.3 ±0.2 


3.2 ± 0.3 


5.1 ± 0.3 


3.9 ±0.4 


C 


1.9 


2.7 ±0.2 


1.8 ±0.3 


3.2 ± 0.3 


2.2 ±0.4 


D 


2.0 


3.6 ±0.2 


2.3 ±0.3 


4.3 ±0.3 


2.7 ±0.4 


E 


3.2 


3.4 ±0.4 


1.7 ±0.5 


4.6 ±0.5 


2.3 ±0.7 


F 


1.1 


5.7 ±0.5 


< 2.2 


6.4 ±0.6 


< 2.4 


G 


1.3 


3.6 ±0.5 


< 2.1 


4.1 ± 0.6 


< 2.3 


H 


1.5 


3.1 ±0.4 


< 1.4 


3.6 ±0.4 


< 1.7 


I 


2.8 


2.8 ±0.3 


1.6 ±0.4 


3.7 ±0.4 


2.1 ±0.6 


J 


2.0 


1.8 ±0.3 


2.0 ±0.4 


2.2 ± 0.3 


2.4 ±0.5 


K 


2.9 


1.5 ±0.3 


< 1.0 


1.9 ±0.3 


< 1.4 


L 


3.0 


1.9 ±0.2 


1.6 ±0.3 


2.5 ± 0.3 


2.1 ±0.4 


M 


1.7 


2.0 ±0.2 


< 1.0 


2.3 ±0.3 


< 1.1 


N 


2.1 


5.1 ±0.2 


1.8 ±0.3 


6.2 ± 0.3 


2.2 ±0.4 





0.5 


5.0 ±0.2 


1.7 ±0.3 


5.2 ±0.2 


1.8 ±0.3 



2 (im cont. c 


S(l)/Bry 


EW (Br 7 ) 


log AT UV 


Age d 


Mass' 


mjy 


ratio 


A 


10 51 s- 1 


Myr 


10 6 M, 


3.24 ±0.03 
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16 ± 1.0 
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29 ± 1.5 
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0.7 
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0.2 
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0.3 
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< 0.5 
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Quoted errors are la and limits are 3<r. 

a Measured in a 3arcsec aperture and corrected for the velocity field (see Section |2.l| ), 
b Corrected for the extinction listed in column 2 using Aj^/Ay = 0.1. 
c Observed flux density, measured in a 3 arcsec aperture. 
d Assuming instantaneous star formation. 



the alignment of Ha, Br7, 10 /im, 2 fim, and CO emission at 
the ILR. 

The radio jet to the North-East extends only 6 arcsec 
from the nucleus and so has not yet impinged on the ring 
(Wilson & Ulvestad 1983), but it is interesting to specu- 
late on what effect it may have if this happens. Models such 
as those proposed by Rees (1989) suggest that once a jet's 
bow shock has passed, clouds will be subjected to an over- 
pressure of ~ M 2 where M is the Mach number. From the 
half-apex angle of the cone which they interpreted as the 
bow wave generated by the jet's motion, Wilson & Ulvestad 
(1983) estimated that M > 3, sufficient to trigger further 
star formation. 

The ionisation cone from the AGN may have an effect 
on the environment to the North-East of the nucleus since it 
extends up to and beyond the ring, and may be assisting the 
star formation by heating the ISM so that the overpressure 
helps to crush the clouds. Similar mechanisms have been put 
forward for jet-induced star formation where the hot plasma 
behind the shock front overpressures the cloud triggering 
star formation (eg van Breugel et al. 1985; Best, Longair 
& Rottgering 1997). We consider whether this mechanism 
may apply here. When clouds first drift in from the galactic 
disc, they will probably have temperatures T ~ 10 K and 
average densities n c ~ 10 3 cm -3 ; after some time in the en- 
vironment of the inner galaxy, they will heat up and perhaps 
coalesce to form bigger, denser clouds with T ~ 100 K and 
n c ~ 10 4 cm" 3 (Goldsmith 1987; Genzel 1991). These clouds 
will have pressures P/k of 10 4 Kcm -3 and 10 6 Kcm -3 re- 
spectively. Once star formation has commenced H n regions 
can provide the overpressure to propagate star formation 
throughout the remaining GMCs. At least some of the gas 
in the ionisation cone photoionised by the Seyfert nucleus 
is in a high excitation state. This has been modelled by 



a number of authors (Evans & Dopita 1986; Bergeron et 
al. 1989; Nazarova 1994) to have T = 3-5 x 10 4 K and 
n c = 0.1-1 cm -3 . The pressure due to this phase is then 
P/k ~ 10 4 Kcm -3 , clearly not high enough to affect the 
GMCs in the current epoch. On the other hand, the soft 
X-ray morphology, which is aligned toward the North-East, 
has been attributed to a nucleus-driven wind with a pres- 
sure of 2-8 x 10 6 K cm" 3 (Wilson et al. 1992). These authors 
suggest this could confine the narrow-line clouds, and would 
also be sufficient to provide some overpressure on clouds in 
parts of the ring, perhaps providing an external stimulous 
for star formation as required in Section W 



7 THE ROLE OF THE NUCLEUS 

Bland-Hawthorn et al. (1997) have proposed a model in 
which the nuclear continuum luminosity incorporates a big 
blue bump with a turnover at ei = 30 eV 

L e (photons s _1 eV -1 ) = k\E~ 2 ^ z exp [— e/ei] 

as well as a power-law with slope a = 1.9 turning over at 
ei = 100 keV 

L e (photons s _1 eV -1 ) = fee - " exp [— e/e2\ 

(fei and k2 are constants) which add together to make a more 
realistic quasar spectrum of the form shown in Fig. ^. This 
luminosity heats grains at the circumnuclear ring which are 
lying within the ionisation cones. The small grains which can 
easily be heated to high temperatures would be responsible 
for the enhanced mid-infrared continuum while the larger 
grains would re-radiate at longer wavelengths. Their model 
predicts continuum fluxes over a range of wavelengths from 
2-100 /im, and was envisaged in order to explain the 10 fim 
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Figure 3. Contour plot of 1-0 S(l) in NGC 1068. Contour inter- 
vals arc as for Fig. 1, and are not corrected for the velocity field. 
The parallel lines mark the extent of the bar at at 0.5 mjy and 
position angle 46° . The arrows delineate the extent of the ionisa- 
tion cone as described by Unger et al. (1992), at PA 11° and PA 
51°. The cone extends beyond the ring. 



Figure 2. Spectra of the two models of the nuclear luminosity, 
scaled to the estimated ionising flux of Q = 8 X 10 54 sec -1 . That 
used by Storchi-Bergmann et al. (1992) is a simple power-law, 
while that from Bland-Hawthorn et al. (1997) has a more realistic 
quasar form, incorporating a power-law with a turnover as well 
as a 'big blue bump' component. 



emission in the ring near the ends of the bar (Telesco & 
Decher 1988). 

In the previous sections we considered only the role of 
star formation for the emission processes in the ring. But the 
continuum-residual image (Fig ^) shows that the strongest 
2/um emission not associated with the bar coincides almost 
exactly with the 10 fim emission, and hence may be excited 
by the same mechanism. We are unable to impose any useful 
constraints on such a model as the predicted total 2 /im flux 
is only 1 mjy, considerably less than that observed at the 
same positions. 



7.1 Circumnuclear H2 Emission 

Fig. ^ shows the alignment of the bar (parallel lines) and 
the ionisation cone (diverging arrows) with respect to the 
observed 1-0 S(l) emission, although direct interpretation of 
the figure should be approached cautiously due to the cali- 
bration difficulties. The dereddened fluxes in Table show 
that the intensity of the Br7 flux is similar at both ends of 
the bar. On the other hand, there is a significant difference 
in the 1-0 S(l) fluxes: that at the North end of the bar is 
almost twice that at the South end. This is reflected in the 
S(l)/Bry ratio which is 0.7-1.0 to the North and < 0.5 to 



the South. Note that the high S(l)/Br7 ratios observed in 
two other regions (J and L) are due to weak Br7 rather than 
strong 1-0 S(l); In every region except A-D the 1-0 S(l) flux 
is less than 2.5 x 10~ 17 Wm" 2 , while in regions A and B it 
reaches 3.9 x 10" 17 Wm" 2 . We propose that there may be 
an additional excitation source for the H2 in these regions. 

Puxley, Hawarden & Mountain (1990) proposed the use 
of the S(l)/Br7 ratio as a diagnostic for models of star 
forming regions, assuming that the 1-0 S(l) was excited by 
fluorescence. Their model B, a compact stellar cluster sur- 
rounded by a single large Hn region, predicts a ratio of 0.1- 
0.4 for electron densities > 100 cm" 3 and photodissociation 
region gas densities 10 3 -10 4 cm" 3 and A?uv = 10 51 -10 52 s" 1 . 
Their model A, involving less compact clusters in which each 
star has its own H 11 region, produces higher ratios for the 
same gas densities. However, there is considerable evidence 
that stars form in compact clusters (eg Meurer et al. 1995), 
which would tend to result in a low S(l) /Br7 ratio for purely 
fluorescent H2 emission, so that for the ratios observed here 
fluorescence is unlikely to be the only source of 1-0 S(l) flux. 

Supernova remnants, appearing in a starburst after 
3 Myr, will produce some emission by shocking the ISM; 
also, where there is such a high density of clouds, shock ex- 
citation directly via cloud collisions may also be an impor- 
tant contributor. Together these may explain the observed 
tendency towards S(l)/Br7 ~ 0.5 . However, all the pro- 
cesses mentioned so far will follow either the star cluster 
or cloud distributions which are fairly uniform around the 
ring. To account for the large difference in 1-0 S(l) flux at 
opposite ends of the bar, we propose that the extra emission 
at the North end, 1-2 x 10~ 18 Wm" 2 in the 3 arcsec aper- 
tures, is due to a process unrelated to star formation and 
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cloud collisions. One possibility is fluorescence at the edges 
of the molecular clouds by UV radiation from the nucleus. 
That the nucleus has an important influence on this region 
is unequivocally demonstrated by the [O in] ionisation cone 
(Unger et al. 1992). Bergeron et al. (1989) showed that the 
emission line spectrum in this region could be modelled in 
terms of photoionisation by soft X-rays from the nucleus. 
We argue that the UV radiation from the nucleus may also 
be affecting the line excitation within the molecular clouds 
in this direction. 

In order to quantify this effect, we convert the incident 
UV flux into an emitted 1-0 S(l) intensity, using the detailed 
models of H2 excitation from Black & van Dishoeck (1987). 
These require us to estimate the 912-1100 A flux, which we 
can do using the same prescription that Storchi-Bergmann, 
Mulchaey, & Wilson (1992) used to calculate the luminos- 
ity of the torus in NGC 1068. The starting point is the ob- 
served characteristics of the ionisation cone, which has a 
direct line of sight to the nucleus. Then by representing the 
nuclear continuum emission as a power-law L v = Au~ a , we 
can derive the coefficient A, and estimate the 912-1100 A 
flux incident at the circumnuclear ring (Fig. Storchi- 
Bergmann et al. calculate the number of ionising photons 
Q from the ionisation parameter U (from [O ill] A5007/H/3 
and [N 11] A6583/Ha) and gas density n (from the ratio of 
the [S 11] lines) measured at several places in the ionisation 
cone at distances r from the nucleus: 

Q = U Anr 2 nc 

The error estimates lead to an uncertainty in nuclear lumi- 
nosity of a factor of 8, but typical values are U = 10~ 3 and 
n = 150 cm 3 at r = 4kpc. The coefficient A is found by 
integrating over all v ^ 1/0, where uo is the frequency corre- 
sponding to the Lyman limit of 912 A, giving A = Qhav^ . 
The spectral index was measured by Kinney at al. (1991) to 
be a = 1.6, similar to the typical Seyfert 1 value of 1.5. The 
912-1100 A luminosity is simply 

Luv = / Av~ a dv — Qhai/Q (vl~ a — i>\~ a )/(l — a) 

J V\ 

where v\ and V2 are the frequencies for 1100 and 912 A re- 
spectively. Similarly, integrating over all v > Vq gives the 
total ionising luminosity as 5 x 10 37 W. 

Sokolowski, Bland-Hawthorn & Cecil (1991) use a simi- 
lar method but proceed in the opposite direction: they begin 
by assuming a luminosity for the nucleus and model the in- 
terstellar medium by fitting the observed line profiles. They 
estimate 10~ 31 ^ U ^ 10~ 2 ' 6 on the ionisation cone and 
10" 4 ' 5 sC U «C 10" 41 away from it for r = 2-4 kpc, finding 
it optimal to set n = 30 cm -3 at 2 kpc, decreasing as r _1 . 
However, their estimate of the nuclear luminosity from X- 
ray observations does not take account of internal extinction 
and ignores the possibility of significant flux from photons 
with E> 15keV. As such it is about 1/4 that of both the 
nuclear infrared luminosity (Telesco et al. 1984) and the esti- 
mated nuclear bolometric luminosity (Pier et al. 1994). If we 
corrected for this, the numerical factor would carry through 
the model and increase their parameter U x n by a similar- 
amount, making it consistent with that of Storchi-Bergmann 
et al. 

Using the values above, we find Luv = 6 x 10 36 W. 
Since there must be rather little extinction in the direction 



of the ionisation cone (we assume none) the UV flux inci- 
dent at the circumnuclear ring, 15.5 arcsec from the nucleus, 
is 4.4 x 10~ 4 Wm~ 2 , nearly 3 orders of magnitude higher 
than the mean background in the solar neighbourhood of 
6 x 10 -7 Wm~ 2 (Black & van Dishoeck 1987). Lastly, if the 
molecular clouds have densities of ~ 10 4 cm -3 and we take 
the covering factor of the clouds to be unity (Planesas et 
al. 1991 find it to be significantly higher than for Galactic 
GMCs) then in a 3 arcsec region the models of Black & van 
Dishoeck give the 1-0 S(l) flux as 4 x 10 -18 Wm~ 2 . 

Two important points to consider are: (1) the factor 8 
uncertainty from Storchi-Bergamnn et al makes little differ- 
ence in practice as in all cases the UV intensity would still 
be more than 100 times the background intensity, and the 
resulting 1-0 S(l) line flux would be reduced by at most a 
factor 2; (2) Sokolowski et al.'s model indicates that the UV 
intensity away from the ionisation cones is 30 times less than 
in them, and is insufficient to cause fluorescent H2 emission. 

We also consider whether this scenario is compatible 
with the quasar-nucleus model of Bland-Hawthorn et al. 
(1997) described at the beginning of the section. The quasar 
spectrum they employed, scaled to the same Q that we cal- 
culated previously as shown in Fig. |^, results in an ionising 
luminosity the same as for the simple power law. But the 
1100-912A luminosity is a factor of 2 less (as the spectrum 
is essentially flat at these energies), although within the un- 
certainty of our previous derivation and sufficient to excite 1- 
S(l) fluorescence at the ring. In order to explain the 10 rim 
emission, Bland-Hawthorn et al. required an EUV (10 eV - 
10 keV) luminosity of 3 x 10 38 W, a factor of 6 more than 
the ionisation and density parameters suggest. However, the 
resulting 1-0 S(l) emission is not significantly enhanced. 

We have demonstrated that in principle a nuclear con- 
tinuum which reproduces the characteristics of the ionisa- 
tion cone can also cause fluorescence at the circumnuclear 
ring, and that this mechanism accounts for both the extra 
flux observed as well as the restricted location. 



8 NEAR-NUCLEAR H 2 EMISSION 

We include a short section here on the nuclear H2 emission 
as our data show, due to considerably lower detection limits 
and wider field, that the emission is more extended than was 
previously realised. Hall et al. (1981) found that the emission 
line spectrum indicated the H2 was excited thermally with 
T cx < 2500 K concluding, as did Oliva & Moorwood (1990), 
that it should be associated with stellar processes. On the 
other hand, Rotaciuc et al. (1991) and Maloney (1997) ar- 
gued that a nonstellar source of ionising radiation such as 
an X-ray emitting AGN, could be responsible for heating 
the clouds. Blietz et al. (1994) suggested that there was a 
close physical connection between the H 2 clouds and radio 
jet, although the nuclear H2 knot could not represent the 
torus proposed by Antonucci & Miller (1985). 

Our observations are not of a sufficiently high resolu- 
tion to shed new light on this debate. However, they do 
show that weak H2 emission extends for 10 arcsec (700 pc) 
along the bar. This is more remarkable because the intensity 
falls away extremely quickly perpendicular to this axis, more 
than an order of magnitude in 1 arcsec. Such a morphology 
is unlikely to arise as a result of collimated emission from 
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Figure 4. Contour plot of 1-0 S(l) in NGC 1068. Contours are 
from 3.6 x 10~ 19 W m 2 arcsec 2 (3.8<r), with levels increasing 
by a factor of 2, to 229 x 10" 19 W m 2 arcsec 2 . These have not 
been corrected for the 250 km s -1 line width. The lowest levels 
have been smoothed to a resolution of 1.7 arcsec, while the highest 
levels are at the full resolution of 1.0 arcsec. The parallel dotted 
lines are the same as those shown in Fig. 3 and delineate the edges 
of the bar as derived from the 3rd contour in the continuum image. 
Those pixels near the nucleus which are above some threshold 
in the on-line image have been blacked out, as these are most 
affected by saturation. 

the nucleus, such as fluorescence due to X-ray irradiation 
of molecular clouds. This mechanism would only be viable 
in the directions along which the high energy photons can 
escape, that is within the ionisation cone. We also consider 
that star forming processes cannot be responsible as the Br7 
flux has a rather different morphology and the S(l)/Br7 ra- 
tio is > 1. Both of these are contrary to expectations of star 
forming regions. 

The alternative we propose is that, since the emitting 
region appears to coincide with the bar, the excitation mech- 
anism is probably connected with it. The bar inside the cir- 
cumnuclear ring is a decoupled secondary and typically the 
co-rotation radius of an inner bar coincides with the ILRs of 
the primary bar. Heifer (1997) measured the inner rotation 
curve and found that if corotation occurs at r — 15 arcsec, 
the pattern speed is 160 kms -1 kpc -1 , and that there may 
be inner Lindbald resonances associated with it in the inner 
hundred parsecs. Consequently, for gas at radii r < 15 arcsec 
inflow can proceed right into the nuclear region. Athanas- 
soula (1992) discusses the response of gas to barred poten- 
tials, showing that shocks occur along the leading edge of 
a bar where the gas flow is close to radial. These models 
do qualitatively describe the inner regions of NGC 1068 in 
which the gas is marginally offset from the bar (to the East of 
the North-East arm and West of the South- West arm, Heifer 
1997), and as such may help in understanding the observed 
H2 at radii r < 15 arcsec which shows weak 1-0 S(l) emission 
along part of the bar, with a morphology entirely consistent 



with the CO map. Even though kinematic studies of the 
inner disc have shown that there may be a general expan- 
sion at velocities 50-100 kms -1 (eg Atherton et al. 1985, 
Planesas et al. 1991), we do not consider that this would 
compromise a scenario in which infall is occuring along the 
bar: clouds form a highly dissipative medium and the loss 
of angular momentum suffered by clouds in the shock would 
be sufficient to overcome any reasonable outward accelerat- 
ing effect, with the relatively high density of clouds in the 
bar resulting in a net infall. Spatially and velocity resolved 
observations of HCN (Tacconi et al. 1994) provide further 
evidence that the bar influences the kinematics of the in- 
ner disc, and that gas is streaming along highly elongated 
orbits. These authors concluded that the H2 in this region 
was shock excited, and that dissipation was leading to influx 
to the nucleus at the rate of a few Mq yr~ and a speed of 
~50kms _1 . 



9 CONCLUSIONS 

We have presented the first near-infrared line images of 
NGC 1068 to include the circumnuclear ring. The fluxes are 
calibrated to take into account the strong variations in the 
velocity field across the inner few kpc, which meant that the 
line flux transmitted by the Fabry-Perot also varied across 
the image. 

We have derived the extinctions in a number of star 
forming knots by measuring the Ha/Br7 recombination line 
ratios. These are typically in the range Av = 1-3 mag. Com- 
parison of the absolute line fluxes and the continuum due to 
young stars (obtained by subtracting the underlying ellipti- 
cal component) to evolutionary synthesis models shows that 
the star formation in the ring occured during a short burst 
within the last 30-40 Myr, and perhaps as recently as 4- 
7Myr ago; and that this was the first such episode of star 
formation in the ring. Although the ionisation cone crosses 
the North-East segment of the ring, it cannot have a sub- 
stantial effect on the star formation. On the other hand, the 
X-ray wind may do so, and once the radio jet reaches the 
ring it is likely that this will also affect it. 

The sizes and masses of the Br7 knots indicate that the 
stars are forming in compact clusters. The 1-0 S(l) fluxes 
and S(l)/Br7 ratios in these regions suggest some compo- 
nent of the 1-0 S(l) emission is due to shock excitation from 
cloud collisions or supernova remnants. However, there is an 
excess flux to the North-East, coincident with the region in 
which the ionisation cone crosses the ring. We have investi- 
gated whether this could be due to fluorescence at the edges 
of molecular clouds by UV radiation from the Seyfert nu- 
cleus which is also producing the high excitation lines, and 
conclude that such a process is possible. 

We have found that the nuclear H2 emission is more 
extended than was previously realised, and that this occurs 
only along the major axis of the bar. It is possible that this 
is associated with shocks induced by the streaming motions 
of clouds along the leading edge of the bar. 
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